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EXCITATION  OF  CURRENT  ON  AN  INFINITE  HORIZONTAL  WIRE 
OVER  THE  EARTH  BY  AN  ARBITRARY  ELECTRIC  DIPOLE  SOURCE 

by 

Ahmed  Hoorfar,  Edward  F.  Kuester  and  David  C.  Chang 


1.  Introduction 


The  propagation  of  discrete  modes  along  an  infinite,  thin,  horizontal 
wire  located  above  the  surface  of  the  earth  has  been  studied  by  many  authors 
([1-3];  see  also  the  bibliography  in  [4]).  Until  recently,  however,  no 
attempt  had  been  made  to  study  the  excitation  of  these  modes  by  anything 
other  than  a delta- function  generator  in  the  wire.  Kuester  and  Chang  [4] 
have  given  formulas  for  the  excitation  coefficients  for  an  arbitrary 
distribution  of  external  current  sources,  while  Wait  [5]  discussed  the 
special  case  of  a vertical  electric  dipole  (VED)  source,  and  also  the  quasi- 
static limit.  Olsen  and  Usta  [6]  also  discuss  this  special  case,  and  study 
the  limiting  case  of  large  distances  away  from  the  dipole  along  the  wire. 

The  somewhat  simpler  situation  of  a perfectly  conducting  earth  is  investigated 
in  [7].  Virtually  no  numerical  results  are  available  for  the  actual  current 
induced  on  the  wire  as  a function  of  dipole  position.  In  fact,  even  for  the 
much  simpler  case  of  a dipole  exciting  a wire  in  free  space,  only  some 
analytical  approximations  for  the  current  at  large  distances  from  the  dipole 
seem  to  be  available  [8,9]. 

The  dipole  excitation  problem  is  of  interest  for  a number  of  practical 


applications.  One  is  the  problem  of  interference  on  open  wire  conanuni cat ions 
from  external  radiators.  Another  might  be  the  distortion  caused  by  a long 
wire  such  as  a power  line  on  the  fields  of  a small  transmitter  used  in  a 


used  in  a direction-finding  scheme.  The  immediate  motivation  for  the 
present  investigation,  however,  is  the  problem  of  mode  conversion  caused 
by  the  proximity  of  a scattering  object  to  a line  propagating  a discrete 
mode.  It  is  possible  that  such  a mechanism  might  be  useful  in  an  intruder 
detection  system  employing  a long  cable  located  in,  on,  or  above  the 
ground  at  the  periphery  of  the  region  to  be  guarded.  An  intruder  disturbing 
the  fields  of  the  mode  normally  propagating  along  the  wire  would  cause  a 
change  in  response  at  a monitoring  station  at  the  end  of  the  wire.  A 
number  of  potential  instrumentation  schemes  have  been  proposed  [10 -l?|. 

In  the  first  part  of  this  report,  based  upon  the  formulation  given 
in  (4),  an  expression  for  the  current  induced  on  a horizontal  wire  above 
the  ground  by  an  arbitrarily  oriented  Hertzian  electric  dipole  is  presented. 
The  total  current  is  the  sum  of  discrete  modal  contributions  as  well  ns 
a number  of  branch  cut  integrals  corresponding  to  radiation  spectra. 

Those  currents  are  evaluated  numerically  for  a variety  of  dipole  orienta- 
tions, and  locations  with  respect  to  the  wire.  It  is  found  that  certain 
configurations  can  excite  the  recently  discovered  earth-attached  mode  [13] 
with  significant  amplitude  compared  to  either  the  ordinary  transmission 
line  (or  structure  at t ached) mode,  or  the  radiation  modes. 

In  the  second  part,  the  scattering  of  a mode  incident  onto  a dielectric 
obstacle  in  proximity  to  the  wire  is  studied.  By  using  a Rayleigh  approxi- 
mation [14]  for  the  induced  polarization  currents  in  the  obstacle,  it-, 
effect  can  bo  replaced  by  that  of  an  electric  dipole  of  appropriate  intensity 
and  orientation.  The  conversion  of  current  into  other  modes  (discrete  as 
well  ns  radiation  modes)  is  examined,  and  the  feasibility  of  this  mechanism 
for  the  detection  of  human  or  animal  intruders  is  discussed. 
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2,  Excitation  by  a Dipole  Source:  Formulation 


The  configuration  to  be  analyzed  is  depicted  in  Fig.  1.  An  arbitrarily 
oriented  electric  Hertzian  dipole  is  located  in  the  air  at  a point  P CXQ »yQ> ZQ) 
as  shown  (xo>0),  above  a conducting  earth  (x  < 0)  with  electrical  parameters 
e,  y and  o . The  dipole  excites  an  infinite  horizontal  wire  of  radius  a 
located  at  a height  x=h  above  the  earth’s  surface,  and  at  y=0.  The 
earth  is  assumed  to  be  homogeneous,  and  a complex  refractive  index  n is 
defined  for  it  as 


n = (e  + ) 

v r (jje 


o 

with  Im(n)  > 0,  appropriate  to  an  exp(-iut)  time  dependence,  which  will 
be  suppressed  hereafter.  It  is  further  assumed  that  a « h and  k^a  « 1, 
so  that  the  thin  wire  approximation,  which  postulates  a uniform  current 
distribution  about  the  circumference  of  the  wire  surface,  can  be  invoked. 

The  integral  equation  for  the  total  current  induced  on  the  wire  by 
a vertical  electric  dipole  source  has  been  derived  by  Wait  [5]  and  by 
Olsen  and  Usta  [6].  A more  general  formulation  valid  for  an  arbitrary  source 
configuration  is  given  by  Kuester  and  Chang  [4].  According  to  their  results, 
the  induced  current  on  the  wire  is  given  as  a Fourier  integral 

r - ikn“2 

1(2)  = KcOe  0 da  (1) 


- CO 


where  kQ  = w/uo£0  is  the  wavenumber  in  air.  Using  a reciprocity 
argument,  the  transform  function  1(a)  is  obtained  as  an  integration  of 
the  product  of  the  source  current  distribution  with  the  fields  produced 
by  a current  on  the  wire.  If  only  external  electric  currents  are  present. 
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In  (2),  Ew  is  given  more  specifically  by 


-ik  a(z-z') 


k t -ik  a(.z-z  j 

Ov-»  ■ i?  I &*’>•  d(I-I') 


(3) 


where  Ew(x;x')  is  the  field  of  a ring  of  axially-directed  electric 
o 

dipoles  of  radius  a,  located  at  a point  x in  the  absence  of  the  wire. 

£W(xt;  -a)  is  thus  the  electric  field  produced  by  an  axial  current  of  the 

form  cxp(-ik  az')  on  the  wire.  J , on  the  other  hand,  represents  the 
*o  c 

external  electric  current  source:  the  arbitrarily-oriented  dipole  in  the 
present  context  which  is  given  by 


J®  (x)  * aepfi(x-xo)6(v-yo)6(z-zo) 

In  all  the  above,  xt  = (x,y)  denotes  the  part  of  x transverse  to  z, 
From  the  geometry  of  Fig.  2,  we  have: 


(4) 


a = a cos0  ♦ a sin6  cos$  ♦ a sin6  sin$ 
ex  y z 


(5) 


while  p = IdJl  is  the  dipole  moment,  and  nQ  = (uQ/e0)  is  the  wave 
impedance  of  free  space. 

The  function  $!{a)  is  the  characteristic  modal  expression  for  the 
wire  (4 J : 

hi  (.a)  = J, 

' l V W V V 

’(<>) 


VA°Y2lH01)(AC)  ‘ »o1)(2Hc)J0(AC)1  * J0(AOtP(a)-a2Q(a)]j 


where  A = kQa  and  H = kQh  are  the  normalized  radius  and  height  of  the 
wire  and  JQ  and  are  the  Bessel  function  and  the  Hankel  function  of 
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Figure  2.  Dipole  orientation 


Figure  3.  Deformation  of  contour  in  the  complex  o-plane 
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first  kind  respectively.  P and  Q are  defined  as  Somroerfei  integrals: 


P(a)  - P(H,0;a)  , Q(a)  » Q(H,0;a) 


where 


P(X,Y;a) 


2 * exp(-u^X+iXY) 

in  Uj+  u2 


Q(X,Y;a)  = i- 


2 J exp(-u^X^iXY) 


n uj  ♦ u2 


wherein 


Uj  * (X2-C2)\  u2  * (X2-^  ; Re(u1>u2)>0 
C = (l-o2)*5  , t = (n2-a2)*5  ; Im(c,t  )>0 


Approximate  expressions  for  P and  Q are  given  in  Appendix  A. 
Inserting  (4)  and  (5)  into  (2) , we  obtain 


where 


»(.>  ■ ^ 

n0  0 M(a) 


= ®e'^0  (x0*y0;-a^ 


* ^ox^X0'y0;'a)coS°  + ^oy^X0,y0’"a^sin9coS^  * ^oz(x0’y0;_a*  sin0sin«l) 

iiW 

The  components  of  tQ  are  given  in  Table  1 (see  [4]).  By  superposition, 
we  may  now  study  a dipole  of  arbitrary  orientation  simply  from  a con- 
sideration of  dipoles  along  the  x,y,  and  z directions  separately. 
Finally,  inserting  (11)  into  (l),  we  have 


9 


ik0a(z-z0) 
e da 


(13) 


Following  [4,15],  if  z-Zq>0,  we  deform  the  contour  of  integration 

from  the  real  axis  of  the  a-plane  upwards  around  three  branch  cuts  eman- 

2 % 

ating  from  a * 1,  a =n/(n  +1)  , and  n,  respectively,  as  shown  in  Fig.  3. 
o 

In  this  process,  residues  at  the  poles  of  the  integrand  (zeroes  of  Jl(a)) 
are  captured.  Assuming  the  earth  to  be  sufficiently  lossy,  the  integral 
around  the  branch  cut  at  a «=  n can  be  neglected  (15],  as  can  the  residue 
at  any  pole  which  has  a large  imaginary  part  of  a.  In  previous  numerical 
studies  [4,  13],  two  poles  with  small  imaginary  part  have  been  found,  so 
that  under  these  conditions,  the  total  current  splits  into  four  contrib- 
utions: 


KZ)  = Ipl(z)  ♦ Ip2(z)  * IB1(Z)  ♦ lB2(z) 


(14) 


where  I , and  I _ are  residue  terms  from  the  poles  at  a , and  a 

pi  p2  K pi  p2 

corresponding  to  discrete  propagating  modes,  referred  to  as  the  structure- 
attached  and  surface-attached  modes.  We  note  that  as  earth  becomes  more 


conducting  and/or  operating  frequency  becomes  lower  the  structure-attached 
mode  is  commonly  known  as  the  transmission-line  mode  because  it  approaches 
the  ideal  TEM-mode  in  the  limit  of  Jn)-*”,  while  Igj  and  are  the  branch 

cut  integrals  around  a = 1 and  a = a , respectively. 

D 

The  two  residue  contributions  have  the  form 


1 

P 


8pni  ^ 

Vo  ft'  (ap) 


,ikoap(z'V 


(IS) 
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radiated  into  the  air--the  "sky  wave"--can  be  expressed  as  an  integral  along 
the  real  axis  from  0 to  1,  in  addition  to  one  along  the  positive 
imaginary  axis: 


Igi  (z)  = ~l~~ 
81  n0k0 


1 ®U°$eo(a)  'Va)^e  (a)  ik0a(z"zo) 

e da 


-l 


- e da* 

?l0(ia' )«,(!«') 


C16) 


where  the  subscripts  "0"  and  "a"  denote  the  argument  of  (. 

The  second  branch  cut  integral,  which  corresponds  to  fields  radiated 
at  various  directions  along  the  ground,  can  be  expressed  as  an  integral 
with  respect  to  a real  variable  by  writing 


(17) 


where  Cg  is  positive  along  the  lower  side  of  the  branch  cut.  Changing 


the  integration  variable  to  c , we  have 

D 


^ 


Va)i!.o(a)-!,o(a)*e,<a)  ‘k0“(z-!0,,;Bd‘B 

~ ‘ 1 


2 2 1* 

where  a * (ag-tg)  and  lm(a)>0.  The  subscripts  0 and  a now  refer 
to  argCD(i.e.,  the  functions  are  evaluated  at  +cD  or  -£„)• 

D D D 

3.  Current  induced  on  the  wire  by  a VKD 

In  this  section  we  present  numerical  results  for  the  current  induced 
by  a vertical  dipole  in  the  xy-plane,  at  x^  =d,  y^  = b.  In  this  case, 


11 


the  induced  current  is  given  by  (13),  with 

Va)  = Kxid’h'~a)  (19) 

and  £*x  is  given  in  Table  1.  Accurate  computational  formulas  for  P,  Q, 

and  their  derivatives  needed  in  ^ (a)  have  been  derived  elsewhere  [161 

e 

and  are  quoted  in  Appendix  A. 

In  the  first  set  of  results  presented  here,  the  wire  height  in  terms  of 
the  freespace  wavelength  X is  h = 0.24X,  the  wire  radius  a = 0.007A, 
and  the  refractive  index  of  the  earth  is  taken  as  n = 5.3  + i0.45.  This 
situation  corresponds  to  a rather  poorly  conducting  earth,  with  e/to’  29 
and  a - 27  millimho/m  at  a frequency  of  100MHz.  The  dipole  strength 
was  taken  in  all  cases  to  be  Id£  = .005  amp-m.  Fig.  4 shows  the  results 
of  excitation  by  a VED  located  on  the  earth's  surface  below  the  wire  (d=0, 
b=0)  as  a function  of  distance  along  the  wire.  Below  about  z/A=  0.1,  the 
integral  (16)  for  1^^,  converges  very  slowly,  and  the  numerical  integration 
is  quite  time-consuming.  It  was,  however,  verified  that  the  total  current 
smoothly  approaches  zero  as  z-*-0,  as  required  by  expression  (13),  since 

is  an  odd  function  of  a.  We  observe  that  the  primary  radiation  current 
Igj  decays  quickly  with  z,  while  1^  drops  off  more  slowly.  Moreover, 
the  pole  term  I (which  corresponds  to  a "structure-attached"  or 
"transmission-line"  mode)  is  initially  the  largest  term,  but  decays  more 
rapidly  than  I (corresponding  to  the  "surface-attached"  or  "earth- 
attached"  mode).  This  latter  mode,  on  the  other  hand,  is  excited  much 
less  efficiently  than  is  the  former,  as  might  be  expected  since  its  fields 
are  not  as  concentrated  underneath  the  wire. 


Fig.  5 shows  the  total  current  and  its  four  components  for  a fixed 


r3  amperes 
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value  of  z as  the  height  of  the  dipole  is  varied.  At  d = h,  the  dipole 
touches  the  wire  and  consequently  the  current  becomes  large  in  this  vicinity. 
As  d increases,  the  first  branch  cut  current  1 decays  quickly  after  a 
few  oscillations,  while  I ^ falls  off  more  slowly.  For  sufficiently  large 
heights,  in  fact,  and  I the  current  of  the  earth-attached  mode, 

are  the  largest  in  magnitude,  but  substantially  cancel  each  other  as  can 
be  seen  from  the  much  smaller  magnitude  of  the  total  current.  For  d 
smaller  than  h,  it  is  clear  that  the  major  part  of  the  current  is  that  of 
the  transmission-line  mode,  which  is  in  agreement  with  the  results  of 
the  previous  figure.  Ne  note  that  the  propagation  constant  a of  the 
transmission-line  mode  for  this  case  is  0.9904b  ♦ iO. 01562,  while  that  of 
the  earth-attached  mode  is  0.99245  + iO. 00239.  The  branch  point 

a is  located  at  0.98301  + iO. 00283. 

B 

In  Fig.  6,  a similar  plot  is  shown  for  a somewhat  larger  value  of 
|n|  and  its  loss  tangent,  n =7.43+  i6.73.  This  case  corresponds  to  a 
lower  frequency  (1MHz,  say),  e/Cq  = 9.9,  and  a rather  smaller  conductivity 
of  o = 5.6  millimho/m.  In  this  case  the  transmission-line  mode  propagation 
constant  = 1.0071  + i 0.0113,  while  that  of  the  earth-attached  mode 

a = 0.9984  + i 0.00371.  The  branch  point  a„  is  now  situated  at  0.99947  + 

p J 1 R 

iO. 00494.  By  comparison  with  Fig.  5,  it  can  be  seen  that  for  this  value  of 
n the  earth-attached  mode  is  excited  somewhat  less  efficiently  when  the 
dipole  is  located  above  the  wire,  indicating  a greater  spread  of  the 
fields  of  this  mode  in  the  region  above  the  wire  for  the  smaller  value  of 
n.  This  slower  decay  rate  also  appears  in  the  slope  of  the  second  branch 


current  I and  can  be  attributed  to  the  decay  rate  of  cxp(-u.  X)  for 


the  associated  fields,  where  u‘ 


2 1 

(n  +1)  . For  Fig.  5 this  value  is 


log  (HI  IxlO6) 


1(> 


Ujp  = 0.00023  ♦ il2.2,  while  for  Fig.  6,  the  decay  rate  is  much  faster  since 
= 0.0044  ♦ il.12.  In  both  cases,  the  wire  has  only  a very  localized 
effect  on  I ,,  while  the  parameters  of  the  earth  play  a major  role  when  the 
dipole  is  not  close  to  the  wire. 

Figs.  7 and  8 examine  the  case  when  the  dipole  is  located  in  the 
xy-plane  at  the  earth's  surface  not  directly  under  the  wire  (x=0,  y - b ) . 

For  a fixed  observation  point  z/A  - 0.5,  Fig.  7 shows  the  effect  of 
varying  b on  the  total  current  and  its  various  components.  For  small  h, 
the  transmission-mode  current  dominates  the  toal  current,  but  falls  off 
rapidly  as  b is  increased.  The  earth-attached  mode  (1^,1  fulls  of  more 
slowly,  as  expected,  but  by  the  time  it  reaches  a magnitude  comparable 
to  that  of  1 ,,  it  is  clear  from  the  value  of  the  total  current  that  a 
good  deal  of  mutual  cancellation  due  to  phase  differences  between  the 
individual  current  contributions  has  taken  place.  In  Fig.  8,  the  currents 
for  a fixed  value  of  b(^l.OA)  are  displayed  as  a function  of  distance 
along  the  wire.  A broad  maximum  in  the  total  current  (due  to  a peak  in 

1 ) occurs  in  the  neighborhood  of  z 0.12A,  but  otherwise  there  is 

o 1 

little  significant  variation  over  a wavelength  or  two. 

In  Figs.  9-11,  the  current  excited  by  a horizontal  electric  dipole 

(HIP)  is  presented,  all  for  the  lower  refractive  index  situation.  Fig.  9 

shows  the  results  for  an  HHP  lying  on  the  earth's  surface  at  b 0.1V 

as  a function  of  distance  along  the  wire.  As  was  the  case  for  Figs.  1 and  $ 

the  t runsmi ss ion- 1 ine  nx>de  current  1 , is  dominant  except  for 

pl  1 

close  to  zero,  when  the  first  branch  current  l becomes  important 

P I 

(compare  Fig.  4).  It  should  be  noted  that  the  total  current  approaches 


zero  as  z ►()  (as  demanded  by  symmetry)  and  that  a broad  maximum  occurs 
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Figure  7 
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b/x 

Current  induced  by  a VCD  on  the  surface  of  the  earth  (x=0.0)  as 
a function  of  the  transverse  displacement  b for  a fixed  distance 
z =0.5X  along  the  wire;  n = 5.3  ♦ i0.45,  h * 0.24\,  a ■ .0071. 


j 
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Figure  8.  Current  induced  by  a VFD  on  the  surface  of  the  earth  at  x/\  - 0.0, 
b/X  * 1.0  as  a function  of  distance  z along  the  wire;  n = 5.3  ♦ 
i0.45,  h = 0. 24\ , a = .007\. 
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around  z = 0.2X.  The  second  branch  current  1^  is  now  much  more 
significant  compared  to  the  earth-attached  mode  current  I ^ than  for 
the  VED.  Fig.  10  shows  the  currents  for  fixed  z and  b as  a function 
of  dipole  height  d.  In  contrast  to  Fig.  5,  the  transmission-mode  current 
I j is  the  larger  part  of  the  toal  current,  and  little  mutual  cancellation 
between  components  of  the  current  occurs,  as  indicated  by  the  magntiude 
of  the  total  current  being  always  larger  than  those  of  its  components. 

As  expected,  a maximum  is  found  when  the  dipole  height  is  equal  to  that 
of  the  wire.  Fig.  11  displays  the  results  for  an  HED  located  on  the  earth's 
surface  as  a function  of  the  transverse  displacement  b.  However,  as  in 
Fig.  7,  substantial  cancellation  between  individual  current  components 
is  indicated  for  large  values  of  b/X.  Note  again  that  all  currents 
vanish  at  b ■»  0,  as  demanded  by  symmetry.  Overall,  the  total  current 
induced  on  the  wire  is  substantially  smaller  than  that  of  a vertical 
dipole.  This  is  particularly  true  when  b/X  is  smaller  than  0.5  wavelength. 


4.  Scattering  by  an  ellipsoidal  dielectric  obstacle 

Of  interest  when  using  the  wire  over  earth  in  a guided  radar  detection 
system  is  the  effect  on  the  wire  currents  of  an  obstacle  in  proximity 
to  the  wire.  Human  or  animal  intruders  are  modelled  as  lossy  dielectric 
scattcrers  (because  of  the  large  water  content  [17])  while  vehicles  or 
other  machinery  must  be  represented  as  perfectly  conducting  bodies.  If 
the  dimensions  of  the  scatterer  are  sufficiently  small,  and  the  incident 
field  sufficiently  uniform  in  the  neighborhood  of  the  scatterer,  a Rayleigh 
or  dipole  approximation  can  be  used  to  determine  the  scattered  fields  [14]. 
Even  though  the  dimensions  appropriate  to  a human  intruder  at  the  frequencies 
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Figure  10.  Current  induced  by  an  MED  at  X = d,  b/X  = 0.1,  z/X  = 0.5  as 
a function  of  height  d;  n = 5.3  + i0.45,  h = 0.24X,  a = .0 
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Current  induced  by  an  III.D  on  the  surface  of  the  earth  at  x/X 
y a b at  r./X  « 0.5,  as  a function  of  b;  n * 5.3  ♦ i0.45. 
h ■«  0.24X,  a «=  . 007X . 


Figure  11. 
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considered  in  this  report  exceed  those  for  which  this  approximation  could 
be  expected  to  be  valid,  it  is  reasonable  to  assume  that  such  an  approx- 
imation can  still  provide  us  with  estimates  of  mode  conversion  due  to  the 
scatterer. 

For  simplicity  we  will  restrict  ourselves  to  consideration  of  an 

ellipsoidal  dielectric  obstacle  with  semi-axes  as  b,  and  c (n>b=c). 

A dielectric  scatterer  can  be  modelled  by  an  electric  dipole  alone,  while 

a metallic  scatterer  requires  a magnetic  dipole  as  well.  If  the  ellipsoid 

is  oriented  along  the  vertical  (x)  axis,  and  an  electric  field  F.  is 

ox 

incident  at  the  center  of  the  obstacle,  then  the  scattered  field  is 
approximately  that  produced  by  an  x-directed  electric  dipole  located  at 
the  center  of  the  ellipsoid  whose  dipole  moment  is  (14J 


where 


Idt 


P c 


4 

'3 


nabc  wP 

x 


P , _!»!!»! — 
X Ij+CN"-!)"1 


(20) 


(21) 


1 ' z!2*5  [*" 


U-!) 


N is  the  refractive  index  of  the  ellipsoid,  and  e is  its  eccentricity 
given  by 


(-’3) 


To  model  a human  body,  we  used  N ~ 12.18  ♦ i 8.27,  corresponding  to 
t - 80  and  o = 0.84  mho/m  at  a frequency  of  7.r>  MHz , with  a/A  ® 0.25 
and  b/A  - c/ A =0.05. 


It 
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Using  this  approximation,  the  mode  conversion  caused  by  an  obstacle 
whose  center  is  located  at  x/A  » 0.25,  y *■  b,  z « 0 when  the  incident 
field  is  that  of  the  structure-attached  (transmission  line)  mode  (a^) 
wus  computed  and  plotted  out  in  Figs.  12  and  13.  The  incident  current 
is  taken  to  be  unity  at  z = 0.  Fig.  12  shows  the  percentage  changes 
in  the  total  current  and  its  components  relative  to  the  amplitude  of  the 
unperturbed  incident  mode  current  in  the  absence  of  the  scatterer  at  an 
observation  point  z/A  ■ 2.0,  while  Fig.  13  displays  the  same  data  for 
z/A  * jo.O,  Although  individual  current  components  can  be  relatively 
quite  large,  it  can  be  seen  that  a good  deal  of  cancellation  occurs 
among  them.  It  should  be  noted  that  surface-attached  mode  1 ,,  at  z - 20. 0A 
has  assumed  a major  importance  because  of  the  smaller  attenuation  of  this 
mode  than  of  1 Additionally,  the  relative  change  in  total  current 
at  b = 0.2A  is  4%  at  20  wavelengths,  but  only  1.3%  at  2 wavelengths. 

Figs.  14  and  15  give  the  corresponding  results  when  the  incident 
mode  is  the  earth-attached  mode  (a  . At  z/A  -2.0  (Fig.  It),  little 
difference  from  Fig.  12  is  seen,  except  for  minor  qualitative  changes 
attributable  to  the  difference  in  field  patterns  of  the  incident  modes. 

The  relative  change  in  total  current  at  b - 0.2X  is  1.01.  At  z = 20. 0A 
(Fig.  15),  these  qualitative  differences  are  exaggerated,  but  it  is 
important  to  note  that  because  of  the  lower  attentuation  rate  of  the 
incident  mode  in  this  case,  the  relative  change  in  total  current  at 
b = 0.2A  is  now  smaller,  about  0.7%. 

We  must  emphasize  again  that  the  obstacle  dimensions  used  here  are 
larger  than  those  for  which  the  Rayleigh  approximation  could  be  con 
side rod  valid;  however,  they  are  not  inordinately  large  and  should 
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Figure  12.  Change  in  current  due  to  scattering  of  transmission  line  mode 
ra  o .99046  + iO. 0156)  from  ellipsoidal  obstacle  located  at 
v P* 

Y - b;  n - 5.2  ♦ i0.45,  h - 0.24*,  a = .007*,  z - 2.0*. 
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Figure  13. 


Change  in  current  due  to  scattering  of  transmission  line  mode 
.99046  ♦ 10.0156)  from  ellipsoidal  obstacle  located  at 


(V 

y ■ b;  n » 5.3  ♦ i0.45,  h * 0.24A,  a ■ .007A  , z 


20. OA. 


Figure  14.  Change  in  current  due  to  scattering  of  earth- attached  mode 

(cjp2  ■ .99245  ♦ iO. 00239)  from  ellipsoidal  obstacle  located  at 

y *=  b;  n n 5.3  ♦ i0.45,  h * 0.24X,  a « .007X.  z = 2. OX. 
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Figure  IS.  Change  in  current  due  to  scattering  of  earth-attached  mode 

(a  * .99245  ♦ iO. 00239)  from  ellipsoidal  obstacle  located  at 

C2bi  n • 5.3  ♦ 10.45,  h - 0.24X,  a « .0071,  z « 20.0X. 
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provide  reasonable  estimates  for  the  scattering  from  obstacles  of  similar 
size. 

5.  Conclusion 

Although  no  dipole  excitation  considered  here  proved  able  to  excite 
the  earth-attached  mode  as  the  dominant  part  of  the  toal  current  on  the 
wire,  this  mode  was  in  many  cases  comparable  to  the  total  current,  and 
suitable  tailoring  of  sources  should  be  possible  such  that  the  amplitude 
of  this  mode  is  enhanced  in  relation  to  the  others. 

Somewhat  paradoxically,  it  was  found  that  an  incident  transmission- 
line mode  provided  a larger  relative  change  in  total  current  at  large 
values  of  z,  at  least  in  part  because  it  attenuates  more  rapidly  than 
does  the  earth-attached  mode.  In  actual  practice,  a localized  generator 
(such  as  a delta-function  voltage  source)  may  be  the  excitation  scheme 
used,  and  both  discrete  modes  as  well  as  radiation  modes  may  well  be 
incident  on  the  obstacle.  Although  the  analysis  of  this  problem  requires 
the  computation  of  radiation  fields  and  not  just  currents  (this  will  be 
taken  up  in  a future  report),  some  predictions  on  the  basis  of  previous 
work  can  be  made.  In  [15]  it  is  shown  that  the  transmission- line  mode 
is  excited  most  efficiently  by  a delta-function  generator  when  the  wire 
height  is  less  than  about  0.15A  above  the  earth.  This  seems  to  indicate 
that,  contrary  to  the  intuition  which  states  that  a large  wire  height 
will  produce  more  field  spread  and  better  detection,  in  fact  a relatively 
small  wire  height  might  be  most  efficient,  in  so  far  as  providing  a 
larger  percentage  of  changes  along  the  wire  at  least  for  a bare  wire. 

It  is  clear  from  Figs.  12-15  that  if  information  about  individual 


Li 


current  components  could  be  obtained,  then  considerably  more  sensitive 
detection  could  take  place.  At  present,  however,  the  authors  are  unaware 
of  such  a selective  detection  scheme,  and  suggest  that  a search  for  one 
would  be  worthwhile. 
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Appendix  A 

Approximate  formulas  for  the  integrals  in  Table  1. 

The  integrals  appearing  in  Table  1 are  P,  Q,  and  the  derivatives  of  Q 
with  respect  to  x and  with  respect  to  y,  where  the  integrals  P and  Q 
are  given  by  (7)  and  (8).  It  is  thus  sufficient  to  obtain  computational  ex- 
pressions for  P and  Q alone,  and  the  remaining  integrals  can  be  obtained 
by  formal  differentiation.  The  desired  expressions  for  P and  Q have 
been  derived  in  [16].  Only  the  approximate  forms  valid  for  2 1 n“  1 (x+d)  “»1 
and  |n|»l  (which  are  essentially  obtained  by  setting  u^  = -in  in  (7) 
and  (8)  [13,16])  will  be  quoted  here: 

rv 

*12 


P<,<'Y;a)  UH,a)W<u>[~ 


*2  J 


2„2 


- Ho‘  «<12» 


(A.l) 


where 


where 


12 

R^2  = [Y^  + (X  H) In  addition 

Q(X,Y;o)  = \ H61)(CR12}  + W(X»YJ«) 
n tm 


W(X,Y;a)  = [Wx(X,Y;a)  + Wo(Y;a)]exp[-i(X  ♦ H)/n] 


(A.  2) 


(A.  3) 


The  function  Wq  is  given  by 


Wo(Y;a)  = ~ exp(icBY)  ♦ S(Y)  ♦ cos (CgY)*^ (a) 


W (a)  = ? ^ * - 

02  CB  (l/n‘-^) 


221--2  3 Un[^  + ( -l2)*]  -&n|cl)  (A. 4) 
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2 2 i 

with  Cg  = (Og-a  ) 1 Im(ag)  i 0,  while  the  explicit  square  roots  in 

WQ2(a)  are  to  have  positive  real  part  and  the  principal  branch  of  the 
logarithm  is  to  be  taken.  The  minus  or  plus  sign  is  to  be  taken  according 
to  whether  arg(C)  = 0 or  » respectively,  on  the  first  branch  cut;  the 
minus  sign  is  to  be  taken  on  both  sides  of  the  second  branch  cut,  and  at 
all  other  points  in  the  first  quadrant  of  the  a-plane.  S is  expressible 
as  an  infinite  series  whose  coefficients  are  defined  recursively: 


* (i;RY),u  -ic  Y URY 

S(Y)  m*  I jjjy  [e  -(-l)me  H.UY) 


IQ(00  = Ho°^  + I Wo^YJH^CCY)  - HjUYIH^Uy)] 

IjUY)  - (cyj-^hJ^ccy)  ♦ ] 

ImUY)  = UYJ^H^UY)  ♦ (■-1)(CY)'2H^I)UY) 
-(m-l)2(CY)'2Iin_2(CY)  ; m > 2 


where  (X)  is  the  Struve  function  of  order  j. 

In  similar  fashion,  the  function  Wx  is  given  by 


(A.S) 


+y2] 


(X ,Y ;a)  *■  -in  £ 
m=0 


(X+H)2m+1  I* 


where  again,  the  coefficients  are  defined  recursively: 


f(o)(6)  = H^ue*) 


f(1)(e)  * jH{n(ce*) 

26*  1 

fImJ(P)  = -G~ 1 [ (m- 1 ) f 1 ^ (0)  ♦ flm"2l(0)l; 


(A. 7) 
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Since  expression  (A. 6)  does  no.t  converge  rapidly  if  Y is  small  compared 


2 2 

to  X ♦ H,  an  alternative  expression  for  W,  useful  when  (X  ♦ H)  > 2Y  , is 


needed : 


W(X,Y:cO  = cos (CRY)W(X,0;a)  - in  £ 
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where 


W(X,0;a)  = exp[-i(X  +H)/n]{-inT(X)  + WQ2(a)  + ) 
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with  I given  by  (A. 5)  and  W -(a)  by  (A. 4).  The  functions  R are 
m oz  m 

defined  recursively 
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Comparison  between  the  formulas  in  this  section  and  the  exact  values  of 
P and  Q has  demonstrated  their  accuracy  for  practical  computations  [16]. 


L 


1 


As  seen  in  Table  1,  it  nay  be  necessary  to  evaluate  ^ or  ^ 
for  certain  dipole  orientations.  Computational  formulas  for  these  quant ites 
can  be  obtained  by  formal  differentiation  of  the  expressions  given  in  this 
Appendix. 
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as  a function  of  height  d for  fixed  distance  z = 0.SX  along 
the  wire;  n 5.3  ♦ i0.45,  h - 0 . 24 X , a * . 007.X. 

Current  induced  by  a VF.D  directly  above  anil  below  the  wire  (h  - 0.0) 
as  a function  of  height  d for  fixed  distance  z - 0.5X  along 
the  wire;  n - 7.43  «•  it*.  73,  h - 0.2X,  a - . 0 1 X . 


Current  induced  by  a VFD  on  the  surface  of  the  earth  (x  ~ 0.0) 

as  a function  of  the  transverse  displacement  b for  a fixed  distance 

z - 0.5\  along  the  wire;  n 5.3  ♦ i0.45,  h - 0.24X,  a = .007X. 

Current  induced  by  a VF.D  on  the  surface  of  the  earth  at  x/\~  0.0, 
b/X  * 1.0  as  a function  of  distance  ?.  along  the  wire;  n •••  5.3  ♦ 
iO.45,  h « 0.24X,  a «•  .007\. 

Current  induced  by  an  1IFD  on  the  surface  of  the  earth  at  x/X  = 0.0, 
b/X  v 0.1  as  a function  of  the  distance  z along  the  wire; 
n =■■  5.3  ♦ i0.45,  h = 0.24X,  a • .007X. 

Current  induced  by  an  HF.D  at  x d,  b/X-  0.1,  r. / X - 0.5  as  a function 
of  height  d;  n - 5. 3 ♦ iO.45,  h = 0.24X,  a - .007X. 

Current  induced  by  an  HF.D  on  the  surface  of  the  earth  at  \/X  = 0.0, 
y ••  b at  z/X  - 0.5,  as  a function  of  b;  n 5.3  ♦ i0.45,  h - l).24\, 
a - .007X. 

Change  in  current  due  to  scattering  of  transmission  line  mode 
(«  - . 0004 (<  * iO.OlSt*)  from  ellipsoidal  obstacle  located  at  b; 

n - 5..’  ♦ iO.45,  h - 0.24\,  a - .007\,  z =■  2.0\. 

Change  in  current  due  to  scattering  of  transmission  line  mode 

(n  , ’ .OOOln  ♦ i0.0l50)  from  ellipsoidal  obstacle  located  at  v • b. 

Pl 

n - 5.3  * iO.45,  It  - 0..’4\.  a - .007X.  z =-  20. 0\. 


(hange  in  current  due  to  scattering  of  earth  attached  mode 
(u  - .00245  ♦ it). 00230)  from  ellipsoidal  obstacle  located  at 


y - b;  n 5.3  « id. 45,  It  - 0.24 \,  a * .007\, 


2.0V. 


t'hange  in  current  due  to  scattering  of  earth  attached  mode 
(u  , - .00245  * it). 00230)  from  ellipsoidal  obstacle  located  at  v - b; 

n - 5.3  * i 0.45,  h 0.24\,  a - .007X.  z - 20. 0\. 


.3*. 


